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how exogenous shocks to the exchange rate affect the yield curves, how bond yields co-move 

in different countries and how the exchange rate is influenced by interest rates, macro- 

economic variables and time-varying bond risk premia. 

Estimating the model with US and German data, we find that time-varying bond risk premia 
account for a significant portion of the variability of the exchange rate: apparently, a currency 

a tends to appreciate when investors expect large capital gains on long-term bonds denominated 

G1 in that currency. A number of other novel empirical findings emerge. 
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1. Introduction 


After the seminal contribution of Ang and Piazzesi (2003), several recent studies have developed no-arbitrage term structure 
models to determine how macroeconomic variables affect bond prices and bond risk premia. While the vast majority of these 
studies have analyzed single countries in isolation, to date very little is known about bond pricing in an international setting and 
how it relates to macroeconomic fundamentals and exchange rate dynamics. 

We propose a two-country no-arbitrage term structure model that can be used to tackle a number of largely unaddressed 
questions about internationally integrated bond markets. For instance, the model allows to assess how exogenous shocks to the 
exchange rate affect the yield curves, how bond yields co-move in different countries and how time-varying bond risk premia 
contribute to exchange rate fluctuations, while also controlling for other macroeconomic variables like inflation and output. 

On the one hand, a number of studies have investigated the joint dynamics of exchange rates, interest rates and other 
macroeconomic variables (e.g.: Eichenbaum & Evans, 1995; Grilli & Roubini, 1995, 1996), without taking bond pricing factors and 
bond risk premia into account. On the other hand, several two-country no-arbitrage term structure models have been proposed to 
analyse the relation between exchange rates, yield curves and bond risk premia (e.g.: Backus, Foresi, & Telmer, 2001; Bansal, 
1997), but in these models all the dynamics are driven by latent variables and macroeconomic variables do not play any role. Our 
paper aims to bridge these two strands of the literature, providing a unified framework to determine how bond-pricing factors in 
two different countries are related to each other and to macroeconomic variables and the exchange rate. 
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Estimating the model with US and German data, some findings emerge that are robust to various specifications of the model and 
choices of the sample period. At short horizons changes in short-term interest rates account for approximately 30% of exchange rate 
fluctuations, while inflation and economic growth have almost no explanatory power. At medium-to-long horizons, the explanatory 
power of inflation and growth becomes higher (up to 25%), while that of interest rates remains about the same. We find that, even after 
accounting for macroeconomic variables and for the impact they have on expectations, uncovered interest rate parity does not hold; 
furthermore, a significant portion of the variability of the exchange rate is accounted for by time-varying bond risk premia (more than 
20% at longer horizons). Conducting an impulse response analysis, we find that increases in bond risk premia trigger a reaction of the 
exchange rate that is similar to the well-known delayed overshooting phenomenon caused by increases in policy rates. Delayed 
overshooting, i.e. persistent currency appreciation after an increase in policy rates, is uncovered by many empirical studies (e.g.: 
Eichenbaum & Evans, 1995; Grilli & Roubini, 1995, 1996) and it is considered one of the puzzles of international finance, as it 
contradicts the theoretical prediction (e.g. Dornbusch, 1976) of an immediate overshooting followed by a subsequent currency 
depreciation. We find that also increases in bond risk premia cause delayed overshooting: a currency tends to persistently appreciate 
when expected excess returns on long-term bonds denominated in that currency rise. As emphasized by Scholl and Uhlig (2006), the 
delayed overshooting puzzle is intimately related to the forward premium puzzle, i.e. the empirical regularity that exchange rate 
fluctuations tend to reinforce rather than attenuate positive return differentials between currencies. According to our estimates, such 
tendency of high yielding currencies to appreciate, found by many researchers with reference to the short-term (and risk-free) 
segment of the bond market (e.g.: Engel, 1996; Fama, 1984), seems to extend also to the long-term segment: when investors expect 
larger capital gains on long-term bonds denominated in a certain currency (in excess of the risk-free rate), that currency tends to 
appreciate. We also find that, after controlling for macroeconomic variables, there are limited spillovers between bond risk premia in 
the two countries, with the result that there is low correlation between bond risk premia in Germany and in the US. Finally, we find that 
exogenous shocks to the exchange rate (those that are not explained by other macroeconomic variables explicitly included in the 
model) have a negligible impact on the yield curves. 

Also Dong (2006) and Chabi-Yo and Yang (2007) have recently studied the behavior of internationally integrated bond markets in 
a no-arbitrage framework with macroeconomic variables. We adopt a modelling strategy which is substantially different from theirs. 
While they take the domestic and the foreign pricing kernel as exogenously given and derive implied currency depreciation 
endogenously, our model features an exogenous process for currency depreciation and an endogenous foreign pricing kernel. This 
approach has several advantages. First, it overcomes the well-known mismatch between model-implied and actual exchange rate 
found when the depreciation process is derived endogenously. This mismatch is commonly attributed to the fact that a predominant 
portion of exchange rate movements is independent of interest rate movements (e.g.: Constantinides, 1992; Leippold & Wu, 2007; 
Lothian & Wu, 2002). The inability to produce a realistic endogenous currency depreciation process is also found by Dong (2006) and 
Chabi-Yo and Yang (2007), notwithstanding the fact that they explicitly consider other macroeconomic factors beyond interest rates: 
both papers find that large variations in the exchange rate remain unexplained by a model with endogenous currency depreciation, 
even after accounting for inflation and output dynamics in the domestic and foreign country. Instead, in our model there is a perfect 
match between the data and the model-implied exchange rate. Furthermore, we allow the exchange rate to be affected by exogenous 
shocks, which can capture important factors not explicitly included in the model, such as current account imbalances (e.g.: Hooper & 
Morton, 1978). Despite having an endogenous foreign pricing kernel, our model is still able to fit very well both the domestic and the 
foreign yield curve. Hence, a more realistic modelling of currency depreciation does not come at the expense of pricing accuracy. 
Another important advantage of our modelling strategy is that it allows to measure the feedback effect of currency depreciation on the 
yield curves: for example, one can estimate impulse-response functions to assess how exogenous currency shocks are transmitted to 
domestic and foreign yield curves and to risk premia. Finally, while in a model with endogenous currency depreciation one has to 
resort to approximate numerical procedures to estimate impulse responses and variance decompositions for the exchange rate (see 
also Diez de los Rios, 2009), our model allows for exact analytical computation of these quantities. 

Although other studies previously recognized that bond returns do not necessarily span returns in the foreign exchange market 
and explicitly introduced exchange rate factors orthogonal to bond market factors (e.g.: Brandt & Santa-Clara, 2002; Graveline, 2006; 
Leippold & Wu, 2007), our study is the first to extend Ang and Piazzesi's (2003) methodology to a two-country setting with exogenous 
currency depreciation. Ang and Piazzesi (2003) have inaugurated a prolific literature which uses modern no-arbitrage pricing models 
to analyze the relation between the yield curve and macroeconomic fundamentals: some examples are Ang, Dong, and Piazzesi 
(2007), Ang, Piazzesi, and Wei (2006), Gallmeyer, Hollifield, and Zin (2005), Hordal, Tristani, and Vestin (2006) and Rudebusch and 
Wu (2004); for a survey, we refer the reader to Diebold, Piazzesi, and Rudebusch (2005). Earlier studies investigating the relation 
between the yield curve and macroeconomic variables, like Fama (1990), Mishkin (1990), Estrella and Mishkin (1995) and Evans and 
Marshall (2007) did not consider no-arbitrage relations among yields and did not model bond pricing. As a consequence, they were 
able to make predictions only about the yields explicitly analyzed (typically no more than three), they did not rule out theoretical 
inconsistencies due to the presence of arbitrage opportunities along the yield curve and, more importantly, they made no predictions 
about risk premia and their dynamics. For these reasons, the more recent studies we mentioned above have proposed to enrich macro- 
finance models with rigorous asset pricing relations, imposing no-arbitrage constraints on bond prices. Our methodological 
contributions to this literature are in the following directions: we propose a general setting for two-country no-arbitrage macro- 
finance models, refining the canonical form proposed by Pericoli and Taboga (2008) for the single-country case; in particular, we add 
further identification restrictions to their canonical form, addressing the lack of identification pointed out by Hamilton and Wu (2010); 
we show how to introduce an exogenous process for currency depreciation in a setting with both observable and unobservable 
variables and we derive a new set of pricing equations which extend the Riccati equations usually found in discrete-time single- 
country models of the term-structure; we develop an OLS-based estimation method that affords considerable simplifications over 
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commonly employed maximum likelihood methods and is better suited to tackle the high dimensionality of two-country models; 
finally, we propose a strategy that allows to identify unobservable bond pricing factors and to map them to easily interpretable 
economic variables. 

The paper is organized as follows: Section 2 introduces the new class of affine two-country models; Section 3 describes the 
estimation method; Section 4 discusses the empirical evidence; Section 5 concludes. 


2. The model 


Our model of the term structure is a discrete-time Gaussian affine model. Let (Q, F, F,, P) be a filtered probability space. 
The dynamics of the vector of state variables X; (a k-dimensional F,-adapted process, k=N) obeys the following stochastic 
difference equation: 


X, = H + pX, + Xe (1) 


where &~ N(0, I.) under P, wis ak x 1 vector and p and 5 are k x k matrices. X can be assumed to be lower triangular without loss 
of generality. The probability measure P is referred to as the ‘historical probability measure’, in order to distinguish it from 
other ’risk neutral’ equivalent pricing measures, to be introduced in what follows. We assume that the first entry of the vector 
X, is the natural logarithm s+; of the exchange rate S, the number of units of local currency you can buy with one unit of foreign 
currency. 

The one-period domestic interest rate r; is assumed to be an affine function of the state variables: 


r, =a+b'x, (2) 


where a is a scalar and b is a kx 1 vector. 
A sufficient condition for the absence of arbitrage opportunities on the bond market is that there exists a risk-neutral measure 
Q, equivalent to P, under which the process X; follows the dynamics: 


X, =—+ PX + Èn: (3) 


where n: ~ N(0, Ig) under Q and such that: 


1. the price at time t of a domestic bond paying one unit of domestic currency at time t+n (denoted by p”) equals!: 
pt = Er [exp(—r,)pr z] (4) 


2. the price at time t of a foreign bond paying one unit of foreign currency at time t+n (denoted by qf?) equals: 
n_ pR Stı n=1 
qe = Ey |exp(—r;) > qt+1 (5) 
t 


The vector ff and the matrix p are in general different from u and p, while ¥ does not change, by equivalence of P and Q. The 
prices of risk, denoted by Ag = }~!(u—ft) and A; = X~! (p—p), provide the link between the risk-neutral distribution Q and the 
historical distribution P: 


dQ 


dP = Fe q 1 / Er [& +1] 


t 


ara] =Ü exp|— (Ao + MX, +j-1)& +3] 


Within this Gaussian framework, bond yields are affine functions of the state variables. For domestic bond yields, the well- 
known pricing formulae (e.g.: Ang & Piazzesi, 2003) are: 


1 
yl = ——In(pl) = An + BX: (6) 


1 EÈ denotes expectation under the probability measure Q, conditional upon the information available at time t. It is a shorthand for E°[.|F;] 
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where yf is the yield at time t of a bond maturing in n periods and A, and B, are coefficients obeying the following simple system of 
Riccati equations, derived from (4): 


A, =a 
B; = 

1 —1 (7) 
An = a fa F (m-1)(A + Bian- BiB, )]| 


Ba = +[b + (n—1)p"B, | 


The pricing formulae for foreign bonds, which we derive in the Appendix A, are: 


1 
z = —_In(qr) = C, + DrX, (8) 


where z/ is the yield at time t of a foreign bond maturing in n periods and C, and D, obey the following system of recursive 
equations: 


Ge af p-3f 22 (9) 
D, =b-pf +f 


1 _ nA 

Ge fa $ n—1)(G4 + EM- ">" ESSE, )| 
1 
n 


[b +f + (n—1)p'E, | (10) 


1 
n—1 = 


fis akx1 vector whose first entry is a 1 and all remaining entries are 0, so that f'X,=s;. 

As clarified by Eqs. (4) and (5), the pricing measure Q can be used to price domestic as well as foreign bonds, with the proviso 
that gains from currency depreciation (S; + 1/S;) must be taken into account to derive the latter. Recent macro-finance studies (e.g.: 
Chabi-Yo & Yang, 2007; Dong, 2006) adopt a different modelling strategy”: the exchange rate is not included in the set of state 
variables X;, the foreign short rate is specified as an affine function of X; and there are two different pricing measures, one for 
domestic bonds (say Q) and one for foreign bonds (say Q’); in this case, assuming market completeness (see e.g. Backus et al., 
2001), exchange rate dynamics are completely determined by Q and Q’. Provided the model is correctly specified, these two 
approaches are completely equivalent: one can either specify the domestic and the foreign pricing kernels and recover the implied 
currency depreciation process, or specify only one of the two pricing kernels together with the currency depreciation process and 
recover the other pricing kernel implicitly (see Graveline, 2006 for a discussion). However, the assumption of correct specification 
is crucial for the two approaches to be truly equivalent: in particular, we prove in the Appendix A that equivalence obtains only if 
the currency depreciation rate (5; +. 1/S¢) is an F-measurable random variable.’ Intuitively, this means that the set of state variables 
included in X; must be sufficiently rich to fully account for the variability of the exchange rate, i.e. the randomness in S; must be 
completely explained by the randomness in X+ This is hardly a realistic assumption when X, does not include the exchange rate as 
one of its components, as shown by the empirical evidence provided by Constantinides (1992), Lothian and Wu (2002), Dong 
(2006), Leippold and Wu (2007) and Chabi-Yo and Yang (2007): although the theoretical exchange rate derived in a model with 
endogenous currency depreciation displays a fair degree of correlation with the actual exchange rate, an important portion of the 
variability of the actual exchange rate remains unexplained. 

We accommodate the presence of an exogenous currency depreciation process, by adopting a modelling strategy popularized 
by Ang and Piazzesi (2003). Their seminal paper proposes a no-arbitrage affine term-structure model in which state variables can 
be both observable and unobservable. A wealth of other papers fruitfully applies their methodology to study the interactions 
between the yield curve and observable macroeconomic variables, for example Rudebusch and Wu (2004), Hördal et al. (2006) 
and Ang et al. (2006). In our model the currency depreciation process is an observable variable and we allow for the presence of 
other observable macroeconomic variables, as well as some unobservables. 

The structure and the parametrization of the model are an adaptation of the canonical form proposed by Pericoli and Taboga 
(2008) for Gaussian homoskedastic affine term-structure models with both observable and unobservable state variables. The 
presentation of the details of the model closely follows Pericoli and Taboga (2008). 


? Diez de los Rios (2009) proposes an approach similar to ours. ` 
3 The filtered probability space (Q, F, F, P) is usually defined implicitly by the law of motion of X;: F, is no larger than the filtration generated by X+, F = o(,0 5 Fi). 
and P is the product measure derived from the transition densities of X;. 
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Suppose that the first k° variables included in the model are observable (remember that the first one is s, the natural logarithm 
of the exchange rate) and the remaining k” = k — k° are unobservable. Collect their values at time t into the k° x 1 vector X? and the 
k"x1 vector X# respectively. Eqs. (1)-(3) can be written as follows: 


Short — rate ï + 
{r, =a + BX? + b" xt 
process 


00 yO OU yu 5 o 
a Et 


Law of motion Xp =p + pXe_y + p™Xi 
(11) 


uU c pu uo yo uu yU uo n0 uu oU 
under P Xp = u" + PPX? + PU Xe + Le + LE 


Law of motion Xp = WD + pMXe_y + POX + VE 


under Q(t apt + exe, + MN + Sn + ait 


where all the matrices are obtained by separating into blocks the matrices in Eqs. (1)-(3). 

The set of identifying restrictions proposed by Pericoli and Taboga (2008) is: 1) °° is lower triangular, 2) }"°=0, 3) 3"“=1, 4) 
b">0, 5) X§=0. To address the potential identification problem pointed out by Hamilton and Wu (2010), we add a further 
restriction: 6) p’” is a block upper triangular matrix satisfying the real Schur property (see the Appendix A for a precise definition). 

The above set of restrictions, defining a canonical form of the model, imposes contemporaneous independence between random 
shocks to the observable variables X? (the exchange rate and other macroeconomic variables) and shocks to the unobservable 
variables Xř related to the shape of the two yield curves. However, as explained in Pericoli and Taboga (2008), this does not imply 
statistical independence of X? and Xř, because the above set of restrictions allows for a lagged response of the exchange rate and 
other observables to changes in the unobservable variables (and vice versa). Also note that the contemporaneous independence 
between shocks to X? and X? does not imply contemporaneous independence between shocks to the exchange rate and shocks to 
the yield curve, because the exchange rate is included in X?, which in turn contributes to determine bond yields via the pricing 
relations (6 ) and (8). Hence, the model allows to reproduce events which contemporaneously affect both the exchange rate and the 
yield curves. 

The above results on identification are easily generalized to the case where the set of state variables includes also some lags of 
the observables; we do not report them here and refer the reader to Pericoli and Taboga (2008). 


3. Estimation method 
In this section we give a detailed account of the estimation procedure we adopted. The less technically inclined reader may 
safely skip this section. 


We use a consistent and exact two-stage estimation procedure, based on ordinary least squares regressions. 
Let 


Y= [yp wave ae yf 
where 7), ..., Ty are the bond maturities used to estimate the model (2N>k"), and let 


A= [Ar, s Ap Cp Cry} 


As we proved in Section 2 and in the Appendix A, in our model yields are affine also in the logarithm of the exchange 
rate: 


Y, =A + BX, (12) 


This property of yields is crucial for the analytical tractability of the estimation method we propose. 
Let T be a k"x1 vector and A a k"x2N matrix (to be specified later). Define 


Then Xf is of the same dimension of X; and follows a vector autoregression: 


Xi =u, + Xia + &.&; (13) 
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where: 


u, = (1-9p07')E + Ou 

p, = Opa! 

3, = a2 
and Q is a (k° +k") x (k° +k") matrix whose first k° rows are the first k° vectors of the Euclidean basis of Re + and the last k" 
rows are equal to the k” rows of AB; & is a (k° +k") x1 vector whose first k° entries are equal to 0 and the last k” equal F + AA. 

Thus, according to Eq. (13) a vector comprising k” linear affine transformation of the yields and the k° observable factors 
follows a vector autoregression; the parameters governing the autoregression are a function of the parameters in the canonical 
form of the model (u, p, >), of the parameters in the affine transformation (T and A) and of the parameters in the pricing equations 
(A and B). (13) is an equivalent representation of the model, as defined by Pericoli and Taboga (2008). In the terminology of Duffie 
and Kan (1996), the representation (13) is a yield-factor model; however, while Duffie and Kan (1996) include only yields in the 
vector of state variables, the vector autoregression (13) includes also other observable variables. 

In order to be able to estimate Eq. (13) one has to fix T and A and assume that the k” linear affine transformations of yields in the 
vector I+ AY; are not subject to mis-pricings or measurement errors. It is standard practice in the term structure literature to 
assume T = 0 and A composed of k” vectors of the Euclidean basis of R2", each selecting only one yield , i.e. k” yields are measured 
without error (Chen & Scott, 1993). We generalize this approach by assuming that the first k” principal components of Y, are priced 
without error. Hence, the rows of A are the eigenvectors associated to the k” largest eigenvalues of the covariance matrix of Y;. We 
make this assumption to maximize the proportion of the variability of yields coming from sources that are not prone to 
measurement or pricing errors. Furthermore, we avoid the somewhat arbitrary choice of which yields are to be considered exactly 
measured and which are not. We also transform the linear combinations thus obtained so that they have zero mean, by setting I = 
— E[AY,]. 

Given the above assumptions and given consistent estimates of IT and A, obtained by calculating the sample principal 
components of Y, , we are able to estimate (13) by ordinary least squares, obtaining consistent estimates of u., p, and $.. 

Now, we are able obtain an equivalent representation X, of X; by rotating and translating X% in such a way that the observable 
variables are left unchanged: 


Xp =, + OX 


where 2, is any (k° +k“) x1 vector whose first k° entries are equal to zero and Q, is any invertible (k° + k") x (k° +k") matrix 
whose first k° rows are the first k° vectors of the Euclidean basis of R® + ™, i.e.: 


Ü, = [ey a Ce Vy we Vw]! 
where: 

e; =[1 0 0 0]! 

e,=[0 1 0 oj 


and vj, ..., Vu are k” vectors such that O» is invertible. 
The equivalent representation X;* has law of motion: 


X” = HU; + Potti F Et 
where 
Ha = (1-2.p,07') E. T OM, 


Pa =O pn)’ 
Da = Q, 5, 


14) 


#, and Q, can be chosen in such a way that u., p.. and X.. satisfy identifying restrictions 1-5 (see previous section). The set of 
restrictions on X., is: 1) ¥2° is lower triangular; 2) 3# =0; 3) 3# =I. 
Since the first k° rows of Q, are the first k° vectors of the Euclidean basis of 
requirement that 2° be lower triangular is trivially satisfied. 
The restrictions >? =0 and }""=1 are satisfied if: 


Re + and 5% is already lower triangular, the 


a a geil, Bie u 
Qo +i =Vj>, i= 1,...,k 
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Since >. is invertible, the restrictions are satisfied with: 


T T -1 
Vi = ep + dy 


Now, if, after performing the above rotations, p*! does not satisfy the real Schur property (restrtiction 6 - see previous section 
and Appendix A), a further rotation is necessary. This is done as follows: 


1. compute the real Schur decomposition of p!": 
ps. = HTH" 


where H is orthogonal and T is block upper triangular with decreasing diagonal entries; 
2. substitute 0+ with RO», where: 


of o0 
R= [o ir] 


The last rotation leaves $» unchanged, because >" = I and H'e}¥ has the same distribution as g#. 

The distribution of any component of €, does not change when you multiply it by — 1, hence you can always change the sign of 
an unobservable component of X;* leaving X.. unchanged, in order to satisfy the restrictions b“ > 0. The restriction Xg*u = 0 can be 
satisfied only by subtracting from the unobservable components of X;* their respective values at t=0. 

The equivalence of X, and Xf and the fact that X;* is an equivalent representation of Xf in canonical form guarantees that 

i; = X+. As a consequence, we can consistently estimate u, p and X in two steps: i) obtain from Eq. (13) consistent estimates of Lt, 
p. and X, by ordinary least squares; ii) recover the coefficients 4.. = 4, p.. =p and X= by continuous transformations which 
preserve consistency (Eq. (14)). 

In the second step of the estimation procedure, we estimate the coefficients A and B in Eq. (12). If Y; denote observed yields and 

E; = Yy, the measurement errors, then: 


Ŷ, = A + BX, + E, (15) 
where, by construction 
AE, = 0 


We propose two alternative estimation procedures: 


1. Ordinary least squares (OLS), i.e. the coefficients A and B in Eq. (15) are estimated running an OLS regression, disregarding the 
no-arbitrage restrictions imposed by the model; 

2. Non-linear least squares (NLS), i.e. the coefficients A and B in Eq. (15) are estimated as non-linear functions of the parameters“f 
and p (so the no-arbitrage restrictions are fully taken into account). 


There are valid arguments in favor of both procedures. The restrictions on A and B imposed by Eqs. (7) and (9) can arguably 
make NLS more efficient than OLS if the pricing kernel is correctly specified (although this is controversial - see Joslin et al., 2011). 
More importantly, i, p and the prices of risk cannot be recovered by OLS, while they can be directly estimated by NLS. On the other 
side of the coin, OLS affords a significant computational advantage, by avoiding the numerically intensive procedures needed to 
solve the non-linear least squares problem; this is especially important given the high dimensionality of the problem at hand. 
Furthermore, OLS avoids the difficulties caused by local minima and flat surfaces encountered by many researchers who estimated 
term-structure models like ours (e.g.: Kim & Orphanides, 2005). Finally, OLS is robust to mis-specification within the exponential 
affine class, in the sense that it produces consistent estimates of A and B also when the pricing kernel is mis-specified, but the true 
pricing kernel gives rise to prices that are affine exponential in the state variables. 

One may argue that setting up a fully-fledged no-arbitrage model is pointless, if one estimates the pricing equations by OLS, 
thereby disregarding the no-arbitrage restrictions. This is not the case, because the no-arbitrage model and the pricing equations 
are indispensable to justify the linear structure in Eq. (15) and to prove that the OLS procedure is indeed legitimate. 


4. Empirical evidence 


In this section we apply the theoretical results outlined in the previous sections to a two-country model where Germany is the 
domestic country and the United States are the foreign country. This section is organized as follows: Subsection 4.1 presents the 
data and briefly discusses some issues related to sample choice and model selection; Subsection 4.2 comments on estimation 
results; Subsection 4.3 examines the ability of the estimated model to accurately reproduce the dynamics of risk premia on bonds 
and currency markets; Subsection 4.4 introduces an equivalent representation of the model that allows for a clear economic 


4 z% has already been estimated in the first step. 
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interpretation of all the state variables; Subsections 4.5 and 4.6 describe the results of variance decompositions and impulse 
response analyses respectively; Subsection 4.7 discusses deviations from uncovered interest rate parity; Subsection 4.8 discusses 
robustness to sample choice; Subsection 4.9 discusses other robustness checks. 


4.1. Data, sample choice and model selection 


We use two datasets of zero coupon rates extracted from US and German government bond yields, recorded at a monthly 
frequency and provided by the Federal Reserve and the Deutsche Bundesbank respectively’: the two yield curves consist of ten 
equally spaced maturities, from 1 to 10 years. The sample goes from January 1973 to September 2007 and the yields are registered 
on the last trading day of each month. We utilize all the ten maturities to estimate our models. 

The macroeconomic variables in our dataset are: an inflation rate and a measure of the output gap for both countries. The 
inflation rate is the twelve-month growth rate of the consumer price index. The output gap is HP-filtered industrial production 
(the smoothing parameter is set to 129,600 as suggested by Ravn & Uhlig, 2002). Both the consumer price indices and the 
industrial production data are taken from Datastream (the codes are BDCONPRCF, USCONPRCE, BDIPTOT.G and USIPTOT.G). We 
also follow Engel and West (2006) and replace the June 1984 outlier in German industrial production with the average of the two 
neighbouring months. 

The dollar/mark exchange rate is the end of month closing middle rate available on Datastream (the dollar/euro exchange rate 
is used to lengthen the time series after December 1998, when the euro became the single currency for the European Monetary 
Union). 

As far as sample selection is concerned, the main discussion in the following sections is based on estimates carried out with a 
sub-sample of the available data, starting from January 1983 and ending in December 1998. This is the same sub-sample used by 
Dong (2006), who provides ample motivation for this choice: according to the evidence presented, among others, by Clarida, Gali, 
and Gertler (1998), it is recommendable to focus on a post-1983 sample in order to avoid shifts in monetary policy regimes both in 
Germany and in the US. The choice of December 1998 as the end date is motivated by the introduction of the euro in Germany in 
January 1999. In Subsection 4.8, where we discuss robustness to sample choice, we summarize the results obtained by enlarging 
the sample to the pre-1983 and post-1998 decades. 

Our preferred model has five observable variables: the logarithm of the exchange rate, German and US inflation and the output 
gap in both countries. The choice of inflation and output gap as the two macroeconomic variables to include is in line with the 
majority of the recent macro-finance literature (e.g.: Ang et al., 2007). The number of unobservable variables is set to four. As we 
explain in more detail in the following sections, four unobservables allow for an excellent bond pricing accuracy and for a neat ex- 
post interpretation of the results in terms of relevant economic variables. No lags of any of the variables are comprised in the state 
vector. In Subsection 4.9 (robustness checks), we consider other specifications where the state vector includes more unobservable 
variables and lags of the observables. 


4.2. Parameter estimates 


Table 1a contains the estimates® of the coefficients of the factor dynamics Eq. (1). The four unobservable variables display a fair 
degree of persistence (the coefficients on their own lags are between 0.865 and 0.975), in line with that of the observable 
macroeconomic variables; the eigenvalues of the autoregression matrix p lie well inside the unit circle, with the lowest modulus 
equal to 0.667 and the highest equal to 0.983, so that the estimated model is covariance stationary. 

The off-diagonal blocks p°" and p"° have many entries that are statistically different from zero (as shown by the high 
t-statistics). As stressed by Pericoli and Taboga (2008), while a number of macro finance studies impose the restrictions 9°” =0 
and p"° =0, these restrictions are not required by the canonical form introduced in Section 2. The overidentifying restrictions 
p*"=0 and p™° =0, together with 5°" = 0 and }“°=0 (which are instead necessary for identification) imply that the observable 
variables are independent from the unobservable and that there are no interactions between macroeconomic variables and other 
factors related to shape of the yield curves (for a discussion of this point, see Rudebusch, Sack, & Swanson, 2006). The relaxation of 
the aforementioned restrictions allows instead for a lagged response of macroeconomic variables to changes in the unobservable 
factors (and vice versa). Recent studies (e.g.: Diebold, Rudebusch, & Aruoba, 2006) find that the hypothesis of independence is 
challenged by formal statistical tests. The evidence presented in Table 1a lends further support to these findings. 

The factor loadings in the pricing equations (15) are estimated both with OLS, exploiting the linearity of the pricing equations 
(see previous section), and with NLS, imposing the no arbitrage restrictions. NLS minimization of the pricing errors is done with 
respect to the prices of risk Ap and Aj, so that these are estimated directly. 

Table 2 reports the estimated standard deviations of the pricing errors of Eq. (15), for both OLS and NLS. Pricing is very accurate 
in both cases, with the standard deviation never exceeding 10 basis points for any given maturity. Notably, foreign bonds are 
priced as accurately as domestic bonds, despite the endogenization of the foreign pricing kernel. 


5 Access to these datasets is granted via the Bank for International Settlements DataBank. 
6 In interpreting results, note that estimation has been carried out with yields and inflation rates expressed in percentage points on an annualised basis. 
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Table 1a 
Dynamics of X,- estimates. 
Equation u Coefficients of lagged variables 
Ex. r. Ger. I. USI. Ger. O. US O. Lt 1 Lt 2 Lt 3 Lt 4 
Exch. r. 0.057 0.916 0.006 0.002 0.001 — 0.006 — 0.003 0.001 — 0.004 0.002 
2.129) (42.96) (2.361) (0.678 (0.581 (—2.256) (— 3.002) (1.766) (— 4.525) 1.891 
Ger. Infl. 0.169 — 0.193 0.935 0.036 0.000 — 0.002 — 0.005 0.005 0.007 0.005 
0.260) (—0.911) (40.21) (1.116 (0.029 (— 0.077) (— 0.476) (0.563) (0.768) 0.477 
US Infl. 0.893 — 0.473 — 0.033 0.907 — 0.021 0.037 — 0.015 0.024 — 0.001 0.0176 
4.115) (— 2.675) (— 1.711) (33.87) — 1.886 (1.680) (— 1.748) (3.562) (— 0.131) 2.23) 
Ger. Outp. — 0.408 0.994 — 0.006 — 0.088 0.702 0.076 0.096 0.027 0.065 0.020 
(— 0.407) (1.214) (— 0.063) (—0.712 (13.54) (0.746) (2.396) (0.874) (1.929) 0.539 
US Outp. 0.793 — 0.564 0.026 —0.057 —0.012 0.904 — 0.029 0.022 — 0.030 0.023 
2.707) (— 2.364) (0.978) (—1.591 —0.778 (30.59) (—2.457) (2.412) (— 3.024) 2.131 
Latent 1 2.327 — 2,339 0.026 — 0.034 —0.018 0.034 0.865 0.048 —0.078 0.042 
2.772) (—3.418) (0.346) (— 0.329 — 0.425 (0.402) (25.94) (1.873) (— 2.768) 1.361 
Latent 2 — 0.691 0.021 — 0.042 0.151 — 0.051 — 0.066 0.000 0.975 0.004 0.002 
(— 0.823) (0.031) (— 0.556) (1.455 = 1.183 (—0.775) - (38.21) (0.153) 0.051 
Latent 3 — 0.008 — 0.323 0.116 0.047 0.082 — 0.106 0.000 0.000 0.932 0.072 
(— 0.009) (— 0.472) (1.546) (0.454 (1.893 (— 1.254) - - (32.74) 2.352 
Latent 4 — 0.403 0.061 0.078 0.026 0.053 0.206 0.000 0.000 — 0.018 0.932 
(— 0.481) (0.089) (1.035) (0.25) (1.232 (2.425) - - (— 0.648) (32.74) 


The table reports coefficient estimates and t-statistics for the vector autoregression Eq. (13). Each row corresponds to a state variable. Each column corresponds to 
each of the lagged state variables on the right hand side of the equations of the autoregression. The column labelled u reports the estimates of the vector of 
intercepts. The sample period is 1983:01 to 1998:12. 


It is interesting to note that the pricing accuracy of NLS is only slightly inferior to that of OLS (6 versus 5 basis points standard 
deviation, averaging across maturities). Furthermore, OLS and NLS fitted yields almost coincide, their differences being on average 
around 3 basis points. This finding could be considered indirect evidence of the fact that the bonds included in our sample are 
priced efficiently, as the deviations of their prices from no-arbitrage levels are economically negligible. This finding also provides 
evidence that in many situations the heavy computational burden imposed by the NLS procedure can be avoided by replacing 
restricted NLS estimates with unrestricted OLS estimates. Of course, the NLS procedure cannot be avoided if the direct object of 
interest are the prices of risk Ag and à4. While we provide comments on the estimated prices of risk in this section, the next 
sections will employ only OLS estimates, so as to illustrate how bond price dynamics and risk premia can be analyzed without 
directly estimating Ap and Aj. 

Table 1b reports the estimated prices of risk. Both Ag and A; have many entries that are statistically different from zero. This fact 
constitutes evidence that the time-variation in risk premia is statistically significant. Furthermore, several off-diagonal entries of 
à are Statistically different from zero, lending support to our unrestricted parametrization (block restrictions on A; such as those 
imposed by Ang & Piazzesi, 2003 would be rejected). 


Table 1b 
Prices of risk - estimates. 
Equation Xo Loadings 
Ex. T. Ger. I. USI. Ger. O. US O. Lt1 Lt 2 Lt 3 Lt 4 

Exch. r. 0.611 — 0.766 0.045 0.031 0.068 — 0.023 — 0.030 0.013 0.023 — 0.046 
(0.021) (—0.036) (0.131 (0.123 (0.421) (— 0.066 (—0.222) (0.016) (0.142) | (—0.100) 
Ger. Infl. 4.188 —9.619 0.717 0.535 — 1.593 0.060 0.371 0.240 0.071 —0.370 
(116.34) (—1.170) (6.940 (7.808 (— 26.92) (0.694 (9.162) (0.926) (1.848) (—2.883) 
US Infl. — 2.615 5.910 — 0.888 0.005 — 0.221 0.547 —0.301 — 0.049 — 0.013 0.389 
— 71.24) (0.376 —3.151 (0.026 (— 1.660) (2.132 (— 2.838) —0.075) (—0.087 1.024 
Ger. Outp. 3.254 — 5.334 0.398 0.133 — 0.091 — 0.787 — 0.330 0.085 0.178 — 0.278 
(759.54) (—3.3924) (18.707 (8.928 (—6.010) (—28.409 (— 33.063) (1.593) (14.349 (— 8.258) 
US Outp. — 3.136 4.217 — 0.072 — 0.222 0.334 — 0.032 0.215 — 0.148 — 0.039 0.1081 
— 82.49) (1.137 — 1.002 (— 5.034 (9.505) (— 0.652 (7.927) —0.962) (—1.204 1.196 
Latent 1 — 3.120 5.316 — 0.515 — 0.102 0.208 0.541 0.027 —0.119 — 0.064 0.279 
(— 679.02) (1.751 (— 10.536 (—3.025 (8.992) (12.053) (1.269) —1.048) (—2.654 4.243 
Latent 2 — 0.546 0.459 — 0.083 0.036 0.027 0.008 0.046 — 0.019 0 0.035 
— 17.26) (0.104 — 0.923 (0.757 (0.644) (0.131 (1.349) — 0.100) (0.001 0.322 
Latent 3 2.293 — 4.636 0.341 0.141 — 0.143 — 0.106 — 0.144 0.068 0.033 — 0.216 
(375.51) (—0.958) (4.891 (2.809 (— 3.786) (— 1.348 (—5.981) (0.416) (0.746 (—2.178 
Latent 4 — 1.652 2.368 —0.097 —0.000 0.407 —0.085 —0.094 —0.079 0.052 0.062 
(— 102.10) (0.555 — 1.389 (— 0.007 (13.98) (— 1.592 (—3.217) — 0.483) (1.843) 0.704 


The table reports the estimates of the prices of risk Ao and A, and their t-statistics. Rows correspond to the state variables whose drift is changed. Columns 
correspond to the state variables that drive the change in the drift. The column Ag reports the time-invariant components of the prices of risk. The sample period is 
1983:01 to 1998:12. 
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Table 2 

Pricing errors. 
Maturity (years) 1 2 3 4 5 6 7 8 9 10 
Panel A - Pricing errors (standard deviations - basis points) 
German bonds - OLS 8 4 5 5 4 3 2 2 3 5 
German bonds - NLS 9 5 6 5 4 3 3 3 5 7 
US bonds - OLS 6 3 4 5 4 4 4 6 5 8 
US bonds - NLS 6 4 5 5 5 5 5 8 6 9 


Panel B - Differences between OLS and NLS estimates (standard deviations - basis points) 
German bonds 5 3 4 3 2 2 2 2 3 5 
US bonds 1 2 2 2 2 4 3 5 3 6 


Panel A reports the standard deviations of the pricing errors of Eq. (15). The standard deviations (expressed in basis points) measure the average distance between 
observed yields and model-implied yields. OLS and NLS refer to the two techniques used to estimate the pricing equations. Each column corresponds to a bond 
maturity. Panel B reports the average distances (measured as standard deviations) between OLS and NLS estimates of model-implied yields. The sample period on 
which the estimates are based is 1983:01 to 1998:12. 


Overall, a qualitative inspection of the estimated prices of risk and of the loadings in the pricing equations does not reveal any 
immediately discernible pattern and the unobservable factors have no straightforward economic interpretation (as, for instance, 
in Ang & Piazzesi, 2003, where the unobservable factors are highly correlated with the level, slope and curvature of the term 
structure). For this reason, we will introduce in the following subsections an equivalent representation of the factor dynamics that 
allows for a clear interpretation of all the state variables of the model. 


4.3. Risk premia 


In this subsection we analyze the ability of our estimated model to reproduce observed patterns of variation in risk premia: we 
perform statistical tests to check whether model-implied expected returns are unbiased predictors of realized returns and we 
assess their in-sample predictive accuracy. We analyze expected returns on domestic bonds and on foreign bonds, considering 
both the perspective of a domestic and of a foreign investor. 

We define risk premia as model-implied expected excess holding-period returns. More precisely, when a domestic bond 
maturing in n periods is bought at time t and held for h<n periods, the risk premium e® " is: 


ef" = E; [In(pt 3h) -In (p;)] + In(p?’) 


The risk premium e% " is the percentage capital gain that an investor expects to obtain by holding an n-maturity zero coupon 


bond for h periods, net of the risk-free interest he could alternatively receive investing in a bond expiring in exactly h periods. 
The risk premium ft " on a foreign bond from the perspective of a foreign investor is: 


fit = E? m(t) inat] + (a) 


while the risk premium g” " on a foreign bond from the perspective of a domestic investor is: 


gl" = E; [In(q 3n) In(qe) + St 4h s,| | In(pr') 


The above definition takes into account the fact that a domestic investor holding a foreign bond for h periods also profits 
(looses) from currency appreciation (depreciation) between t and t+h. 

Realized excess holding-period returns are defined accordingly, substituting expectations of random variables with their 
respective realizations in the above formulae. For example, the realized excess returns re” " on domestic bonds are defined 
as: 


reh” = In(pe 3'r) —In(pr}) + In(pr') 


Table 3 reports the summary statistics of a battery of Mincer and Zarnowitz (1969) regressions of realized excess returns on 
their model-implied expectations: 


ren" =a + Ber" + u, (16) 


We consider both domestic and foreign bonds from the perspective of a foreign and a domestic investor. The holding period is 
one year (h=12 in our monthly model). 
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Table 3 

Regressions of realized on model-implied expected excess returns. 
Bond maturity 1 3 5 10 
German bonds 
R? _ 0.456 0.428 0.359 
a (t-stat) = — 0.915 (— 1.652) — 1.580 (— 1.470) —2.491 (— 1.017) 
B (t-stat) E 1.550 (5.850) 1.568 (5.307) 1.591 (3.710) 
z2 [a=0, B=1] (p-val) _ 4.328 (0.115) 3.723 (0.156) 1.958 (0.376) 
US Bonds 
R? _ 0.355 0.316 0.317 
a (t-stat) = — 0.3817 (— 0.647) — 0.856 (— 0.777) —3.337 (— 1.256) 
B (t-stat) = 1.275 (3.492) 1.350 (3.185) 1.743 (3.512) 
H° [a=0, B=1] (p-val) _ 0.609 (0.738) 0.779 (0.678) 2.360 (0.307) 
US bonds for a German investor 
R? 0.409 0.437 0.444 0.409 
a (t-stat) 0.033 (0.014) — 2.099 (—0.94) — 3.755 (— 1.699) — 4.830 (—2.091) 
B (t-stat) 1.181 (5.259) 1.352 (5.763) 1.419 (5.860) 1.243 (5.447) 
X [a=0, B=1] (p-val) 0.741 (0.690) 2.493 (0.288) 4.564 (0.102) 4.810 (0.090) 


The table reports estimates (t-statistics in parentheses) of Eq. (16), together with R?, and 7” statistics (p-values in parentheses) for the null of unbiasedness. The 
equation is a linear regression of excess bond returns observed during the sample period on a constant and expected excess returns derived with our no-arbitrage 
model. Newey-West standard errors with truncation at 18 lags are used to take into account the serial correlation induced by overlapping returns. The sample 
period is 1983:01 to 1998:12. 


The R? is always higher than 30% and in some cases it is as high as 45%. The estimated Ps are always significantly different from 
zero at all conventional levels of confidence.’ If model-implied expectations are not biased, tests of the joint hypothesis œ= 0, 
B=1 should not be rejected by the data. We run Wald tests of this hypothesis ( y? statistics and p-values are reported in Table 3). 
The hypothesis of unbiasedness is never rejected at 99 and 95% confidence, while it is in only one case rejected at the 90% level. 

The fit of the above regressions is comparable to that found with popular predictive regressions: for instance, Cochrane and 
Piazzesi (2005) also find that their statistical model of excess returns produces R? as high as 45%. Furthermore, the regressions 
provide good evidence that the model produces unbiased estimates of risk premia. As we will explain in the next subsections, 
estimated risk premia play a key role in the analysis and interpretation of model dynamics, hence we believe that analysing the 
properties of regression Eq. (16) is an essential check on the robustness of the model. 


4.4. Equivalent representations and interpretation of the state variables 


In this subsection we introduce an equivalent representation of the factor dynamics (1) that we will exploit in the following 
subsections to carry out variance decompositions and impulse response analyses. While the unobservable variables, as identified 
in the estimated canonical form, are difficult to give a clear economic interpretation, the equivalent representation we now 
introduce allows to map all the state variables to easily understandable economic quantities. We believe this is an important step 
towards a better understanding of the mechanics of a reduced-form no-arbitrage model like ours® and that it helps to provide a 
deeper insight into the economic forces driving bond prices and currency depreciation. 

Note first that the yields on the one-year domestic and foreign bonds are affine functions of the state variables: 


12 
Y = Ay + BioX, 


12 _ T 
Z“ = C2 + DiX: 


As shown in the Appendix A, also the expected excess holding period returns on the ten-year domestic and foreign bonds (the 
latter from the perspective of a foreign investor) are affine in the state variables: 


120,12 _ T 
er = Ji2012 + K120,12Xt 
120,12 _ T 

f = Qy20,12 + Wi20,12X¢ 


where J120, 12 and Q120, 12 are scalars and K120, 12 and W120, 12 are kx 1 vectors whose functional dependence on model parameters 
is reported in the Appendix A. 


7 We use Newey-West standard errors with 18 lags to take into account the serial correlation induced by overlapping returns, as suggested by Cochrane and 
Piazzesi (2005). 

8 To our knowledge, our paper is the first to utilize equivalent representations to achieve economic identification of the latent variables in a reduced-form no- 
arbitrage model. We thank Pietro Veronesi, who encouraged us to think deeply about economic identification of the latent variables. 
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Let T; be a k"x1 vector and A; ak" x (k° +k") matrix defined as follows: 


Ty = [A2 Ci J120,12 Q120:19| 
T 
A, = [Biz Di2 K120,12 Wi20.12| 


Define 


# _ [yoT (12 (12 (120,12 120,12]? 
x = |x; Yt a & fi | 


Then X? is of the same dimension of X, and follows a vector autoregression: 


XP = My + pyXt-1 + Lye (17) 
where: 


He = (I-Ogp05") By + Ogu 


Px = gpg" 

Dy = Q2 
and Oy is a (k° +k") x (k° +k") matrix whose first k° rows are the first k° vectors of the Euclidean basis of RÉ + and the last k" 
rows are equal to the k” rows of 41; Sy is a (k° +k") x1 vector whose first k° entries are equal to 0 and the last k” equal to T4. 

(17) is an equivalent representation of the model, where the state vector X? comprises the 5 initial observable variables, the 

two yields on the one-year zero coupon bonds (the domestic and the foreign one) and the two risk premia® (the expected excess 
holding-period returns on the domestic and the foreign ten-year bonds). Hence, the equivalent representation Eq. (17) allows to 
express the dynamics of the model in terms of easily understandable economic quantities, without resorting to latent factors that 
are not economically identified: as we show in the following subsections, this result allows for a neat interpretation of variance 
decompositions and impulse response analyses. 


4.5. Variance decompositions 


In this subsection we comment on the results of variance decompositions, based on the equivalent representation Eq. (17). 

Table 4 reports the decomposition of the variance of the exchange rate. The proportion of variance explained by inflation and 
output increases with the time horizon, from 6% at a 12-month horizon to 24% at a 120-months horizon. The variance explained by 
shocks to the one-year interest rates is 29% at a 12-month horizon and only slightly decreases as the forecast horizon increases. The 
percentage of variance attributable to risk premia is about 9% at a 12-month horizon and increases up to 24% at longer horizons. 
The proportion of variance of the exchange rate not explained by other variables is at first high (55%), but then decreases to less 
than 25%. In all cases, the predominant role is played by US variables, while shocks to German variables seem to have a fairly 
limited impact on the exchange rate. 

Our findings are in line with those of previous studies which found that, although a significant part of exchange rate 
fluctuations is not explained by macroeconomic variables and interest rates in the short run, over longer horizons the explanatory 
power of fundamentals significantly increases (e.g.: Bekaert & Hodrick, 1992). However, a novel finding emerging from the 
variance decomposition is the fact that a significant portion of exchange rate variability is driven by bond risk premia. As we more 
extensively discuss in the sections on robustness checks, this finding does not appear to be dependent upon any particular 
modelling or identification assumption: in fact, if we estimate a vector autoregression like Eq. (17) that excludes the two variables 
related to bond risk premia, the proportion of variance of the exchange rate left unexplained increases by more than 10% at all time 
horizons. In the following subsections, where we analyse impulse response functions and deviations from uncovered interest rate 
parity, we describe the economic mechanism underlying the influence of risk premia on the exchange rate. 

Table 5 reports the results of decomposing the variance of risk premia at different time-horizons. Results differ depending on 
the country under consideration. On the one hand, about 80% of the variance of the risk premium on the US bond is explained by 
shocks to macroeconomic variables and interest rates, both at short and at medium-to-long horizons. On the other hand, more 
than 70% of the variance of the risk premium on the German bond remains unexplained by macro-factors at short horizons. 
However, also for Germany at longer horizons the explanatory power of the other state variables increases up to about 50%. 
Interestingly, after accounting for macroeconomic shocks, there seem to be very limited spillovers from US to German risk premia 
and vice versa: in each country, the proportion of variance of the domestic risk premium explained by shocks to the foreign risk 
premium is negligible. This is confirmed also by the impulse response analysis in the next subsection. 


? Other economically identified equivalent representations are of course possible. For example, an equivalent representation could be obtained by replacing 
risk premia with term premia. 
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Table 4 

Decomposition of the variance of the exchange rate. 
Proportion of the Forecast horizon 
variance accounted 
for by shocks to: 12 months (%) 36 months (%) 60 months (%) 120 months (%) 
German inflation 1.48 1.19 2.44 3.67 
US inflation 0.66 4.57 6.26 7.38 
German output gap 0.81 4.03 4.61 4.86 
US output gap 2.74 5.28 7.27 7.65 
German interest rate 2.87 4.52 3.70 3.26 
US interest rate 27.17 27.96 25.19 24.96 
German risk premium 2.19 1.28 1.72 1.90 
US risk premium 6.98 18.93 21.85 22.52 
Exchange rate 55.11 32.23 26.97 23.78 


The table lists the contribution of each of the state variables in the vector autoregression Eq. (17) to the forecast variance of the exchange rate, for various forecast 
horizons. The Choleski ordering for the variance decomposition is the same as the order of appearance of the variables in the table (from German inflation to the 
exchange rate). The sample period on which estimates are based is 1983:01 to 1998:12. 


Table 5 
Decomposition of the variance of risk premia. 
Forecast horizon Exchange rate (%) Macro variables (%) Interest rates (%) German risk pr. (%) US risk pr. (%) 
Proportion of variance of the German risk premium explained by shocks to 
12 months 1.08 16.17 7.35 67.88 1.54 
36 months 1.22 30.29 8.45 58.20 1.83 
60 months 1.87 29.88 11.50 54.93 1.82 
120 months 2.92 29.46 13.76 52.02 1.84 


Proportion of variance of the US risk premium explained by shocks to 


12 months 0.12 64.69 19.12 0.50 15.56 
36 months 0.20 52.93 26.41 0.39 20.07 
60 months 0.38 47.59 31,55 0.33 20.16 
120 months 0.90 43.14 35.86 0.36 19.74 


The table lists the contribution of the state variables in the vector autoregression Eq. (17) to the forecast variance of bond risk premia, for various forecast horizons. 
The Choleski ordering for the variance decomposition is the same as the order of appearance of the variables in the table (from the exchange rate to US risk 
premia). The sample period on which estimates are based is 1983:01 to 1998:12. 


4.6. Impulse response functions 


4.6.1. Responses of the exchange rate 

Fig. 1 plots the responses of the logarithm of the exchange rate to one standard deviation shocks!° to output gaps, inflation 
rates, one-year interest rates and risk premia. 

An increase in the German (US) interest rate leads to an appreciation of the Mark (dollar). The maximal response is achieved 
after 13 months for Germany and after 9 months for the US. This is in line with many empirical studies (e.g.: Eichenbaum & Evans, 
1995; Grilli & Roubini, 1995, 1996) that find evidence of persistent currency appreciation after an increase in policy rates. The 
phenomenon is known as delayed overshooting and it is considered one of the puzzles of international finance, as it contradicts the 
theoretical prediction (e.g. Dornbusch, 1976) of an immediate overshooting followed by a subsequent currency depreciation. As 
emphasized by Scholl and Uhlig (2006), the delayed overshooting puzzle is intimately related to the forward premium puzzle we 
discuss in Subsection 4.7. 

An increase in German inflation initially leads to a small depreciation of the Mark, followed, after 6 months, by an appreciation 
of greater magnitude: the overall outcome is a prolonged increase in the value of the Mark with respect to the US dollar. An 
increase in US inflation causes a persistent appreciation of the Dollar. Hence, in both countries, the response of the exchange rate to 
inflation is inconsistent with long-run PPP; it is, however, consistent with the previously mentioned delayed overshooting 
phenomenon: as noted by Clarida (2004) and Dong (2006) central banks tend to increase the interest rate in response to higher 
inflation; if a persistent currency appreciation follows such increase, then the combined effect of these two reactions leads to a 
positive lead-lag relationship between inflation and currency appreciation. 

The response to a positive shock to output is different in the two countries. An increase in German output is followed by an 
appreciation of the Mark, while an increase in US output leads to a depreciation of the dollar. The response of the exchange rate to 
changes in output is notoriously difficult to predict and interpret (e.g.: Harberger, 2003; McFarlane & Zhu, 2006; Miyajima, 2005) 
and the empirical evidence is mixed (Froot & Rogoff, 1994); the response depends on many factors: for example, it may depend on 
whether the output shock is demand or supply driven and whether it has been generated in the tradeables or non-tradeables 
sector; furthermore, it may also depend on how the central bank reacts to output shocks if it is following a Taylor (1993) rule. 


10 We present responses to non-factorized innovations, i.e. for each variable subjected to a shock we plot the effect of a ceteris paribus increase in that variable 
on the other variables. 
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Fig. 1. Responses of the exchange rate to non-factorized one standard deviation shocks to the other state variables. Time (expressed in months) is reported on the 
x-axis of each plot. Thin lines are 90% confidence bands. 
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Fig. 2. Responses of risk premia to non-factorized one standard deviation shocks to the other state variables. Time (expressed in months) is reported on the x-axis 
of each plot. Thin lines are 90% confidence bands. 
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An increase in expected excess holding period returns on German (US) 10-year bonds triggers an appreciation of the Mark 
(dollar); the response is much stronger and more persistent for the United States than for Germany. Hence, the tendency of high 
yielding currencies to appreciate, found by many empirical studies with reference to the short-term (and risk-free) segment of the 
bond market (e.g.: Engel, 1996; Fama, 1984), seems to extend also to the long-term segment: when investors expect larger capital 
gains on Mark (dollar) denominated 10-year bonds (in excess of the risk-free rate), the Mark (dollar) tends to appreciate. As 
discussed also in Subsection 4.7, this evidence adds to the complexity of the delayed overshooting and forward premium puzzles: 
exchange rate fluctuations seem to reinforce rather than attenuate positive return differentials between currencies, not only with 
reference to investments in short-term deposits, but also with reference to longer-term securities. 


4.6.2. Responses of risk premia 

Fig. 2 plots the responses of the risk premia (defined as excess expected holding period returns on the 10-year bonds) to one 
standard deviation shocks to the other state variables. 

A positive shock to German (US) output leads to a reduction of the risk premium on the German (US) bond. The effect is more 
pronounced for Germany than for the US. Countercyclical variation in risk premia is consistent with the predictions of 
intertemporal asset pricing models (Robotti, 2002): in the presence of positive risk aversion (hence decreasing marginal utility of 
consumption), high risk premia are required during unanticipated recessions to induce investors away from (already low) current 
consumption and into risky investments. 

The response to inflation is positive both in Germany and in the US. The response in Germany is very weak and reverts to negative 
after few months, while it is stronger and lasts longer in the US. Hence, investors tend to require higher risk premia when inflation is 
higher; this is consistent with the estimate of a positive price of risk for inflation in Ang and Piazzesi's (2003) Macro Lag Model. 

The risk premium on the German 10-year bond tends to rise when the premium on the analogous US bond increases. Also the 
premium on the US bond has a positive response to increases in the German premium, although to a much lesser extent. The 
correlation between the two is, however, quite low (4%). 

Finally, we do not find any significant response of risk premia to unexpected currency appreciation and to unexpected 
increases in the short-term interest rates. 


4.7. Deviations from uncovered interest rate parity 


The theory of uncovered interest rate parity states that the currency of the country with the higher interest rate tends to 
depreciate, so that the positive yield differential is compensated by an expected currency loss and a risk-neutral investor is 
indifferent between investing in the domestic currency and investing in the foreign currency. There is ample and robust empirical 
evidence (e.g.: Engel, 1996; Fama, 1984) that uncovered interest rate parity does not hold in practice and that higher-yielding 
currencies display a tendency to appreciate rather than depreciate. This phenomenon, widely known as the forward premium 
puzzle, is exploited by a popular trading strategy, known as the carry-trade, whereby one borrows money in low-yielding 
currencies and invests it in high-yielding currencies, hoping to profit both from interest rate differentials and from currency 
appreciation. In this subsection we thoroughly investigate the phenomenon through the lens of our model. 

We need to introduce a fourth definition of risk premium, called foreign exchange risk premium: 


Be" = Et [In(ae sn) —In(ar) + Se + ns] Er [In(pP 3'r) —In(pP) | (18) 


It is the extra-return that a domestic investor expects to gain, over an holding-period of length h, when he substitutes a 
domestic bond maturing in n periods with a foreign bond of equal maturity. Stated differently, it is the expected return from a 
carry-trade strategy whereby one buys a foreign bond and contemporaneously sells shorts a domestic bond of equal maturity. 
While the majority of papers dealing with the forward premium puzzle have limited their analysis of the foreign exchange risk 
premium to the case in which n=h (holding period and bond maturity coincide), some recent papers (Alexius, 2000; Meredith & 
Chinn, 1998) have shown that the empirical evidence is significantly influenced by the choice of maturity. For this reason, we 
choose a more flexible definition of foreign exchange risk premium, where the maturity n is allowed to vary. 

In the notation of our model, the most classical version of the uncovered interest rate parity can be expressed as follows: 


Er [se + n—S:) = In(ar) —In(p7) (19) 
= h(yt —z) (20) 


which is equivalent to: 


gr" =0 


Hence, the main prediction of the uncovered interest rate parity theory is that the foreign exchange risk premium (as 
defined above) be equal to zero, not only on average (i.e. unconditionally), but also in any time period (i.e. conditionally). The 
prediction is a direct consequence of assuming risk-neutral investors. It can be proved that in our model, instead, gia is an affine 
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Table 6 

Affine representation of foreign exchange risk premia. 
Maturities of the bonds used to set up the carry trade 1 3 5 10 
Constant 13.77 11.03 9.66 13.77 
Exchange rate — 49.37 — 47.86 — 47.90 — 49.37 
Domestic inflation =0.30 0.56 0.83 —=0.30 
Foreign inflation 4.93 4.68 4.44 4.93 
Domestic output gap — 1.46 = 1.26 = 126 — 1.46 
Foreign output gap —6.59 —5.80 —5.72 —6.59 
Domestic 1-year interest rate —2.28 — 1.94 — 1.67 —2.28 
Domestic 10-year risk premium 1.89 1.52 1.45 1.89 
Interest rate differential 2.79 3.18 3.45 2.79 
Risk premium differential 2.19 2.11 2.29 3.19 


The table reports the loadings on the set of state variables in the affine representation of foreign exchange risk premia (Eq. (18)). Each column in the table refers to a 
different maturity. The sample period on which estimates are based is 1983:01 to 1998:12. 


function of the state variables and is therefore allowed to vary over time. More specifically, gan is equal to a constant plus a linear 
combination of the following variables'': 1-5) the five observable variables included in X?; 6) the domestic one-year interest rate; 
7) the risk premium on the domestic ten-year bond; 8) the differential between the foreign and the domestic one-year interest 
rate (henceforth, the interest rate differential); 9) the differential between the foreign and the domestic risk premium on the ten- 
year bond (henceforth, the risk premium differential). 

Table 6 reports the coefficients of the linear combination for various maturities (h is fixed to 12 as in the rest of the paper). 
Contrary to the uncovered interest rate parity hypothesis, the foreign exchange risk premium is an increasing function of the 
interest rate differential across all maturities, i.e. a widening of the difference between the one-year foreign interest rate and the 
corresponding domestic rate increases the expected return of the carry trade, irrespective of the maturity of the bonds involved in 
the trade. Note, however, that this is a ceteris paribus statement and does not imply that a positive interest rate differential is 
associated to a positive foreign exchange risk premium, as implied by naive formulations of the forward premium puzzle; in fact, 
the interest rate differential is only one of the variables which concur to determine the foreign exchange risk premium and it may 
well happen that a positive interest rate differential is more than compensated by the negative effect of other variables so as to 
produce a negative foreign exchange risk premium. For example, the coefficients in Table 6 indicate that the stronger the US dollar 
is with respect to the Mark, the lower is the expected profit on the US carry trade (borrowing money in Germany and investing it in 
the US); hence, even when the dollar is the higher yielding currency, the expected return on the carry trade can be negative if the 
USD/DEM exchange rate is sufficiently far above its equilibrium value. 

The foreign exchange risk premium is also an increasing function of the risk premium differential, i.e. the expected return on 
the carry trade increases when the difference between the risk premium on the foreign 10-year bond and the risk premium on the 
corresponding domestic bond widens. This is consistent with the finding reported in Appendix A that higher returns on foreign 
bonds tend to be reinforced by currency appreciation also at longer maturities. 

The dependence of foreign exchange risk premia on other variables is hard to interpret: we report it without further comments. The 
effect of domestic inflation on the foreign exchange risk premium depends on the maturity under consideration, being positive for short 
and long maturities and negative for intermediate. The effect of foreign inflation is unambiguously positive. The effect of the output gap 
(both domestic and foreign) is negative. The effects of the domestic interest rate and risk premium are negative and positive respectively. 


4.8. Extending the sample period 


In this subsection we briefly illustrate the results obtained by estimating the model over longer time spans. We consider two 
prolongments of the 1983-1998 sample period discussed in the previous sections. The first extended sample goes from January 
1974 to December 1998. This is the same period chosen by Inci and Lu (2004), in order to cover the entire floating exchange rate 
period before the introduction of the euro (more precisely, the international exchange rates became floating in August 1973, but 
the authors recommend starting from January 1974 to allow for an adjustment period). The second extended sample goes from 
January 1983 to September 2007 (the month when we last updated our dataset). The exchange rate of individual currencies of the 
European Monetary Union vis à vis the US dollar has been extended beyond 1998 using the dollar/euro exchange rate also by 
Gadea, Montanés, and Reyes (2004), who argue that the extension might augment the power of PPP and unit root tests and that it 
offers a number of interesting insights by reflecting the impact of the early years of the European single currency on the results 
obtained for the pre-euro period. Of course, this approach is prone to several theoretical inconsistencies, but it has the advantage of 
allowing empirical exercises that would otherwise be infeasible. 


11 The proof that the foreign exchange risk premium is an affine function of the variables listed here is very similar to the proof that expected excess holding 
period returns are affine in the state variables. We do not report it to save space. 
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Table 7 

Structural break tests for individual equations. 
Equation Break after 1982 Break after 1998 

Chow B.P. CUSUM Chow B.P. CUSUM 

Exchange rate Yes No Yes Yes 
German inflation No Yes No No 
US inflation Yes Yes Yes No 
German output gap Yes No Yes No 
US output gap Yes Yes Yes No 
German interest rate Yes No No No 
US interest rate Yes No No Yes 
German risk premium Yes No Yes No 
US risk premium No No Yes No 


The table reports the results of test for structural change for the single equations of the vector autoregression (17). The null hypothesis is: no structural change. The 
confidence level is 95%. Yes indicates a rejection, No a failure to reject the hypothesis of stability. The test for a break after 1982 is conducted on the 1974:1998 
sample. The test for a break after 1998 is conducted on the 1983:2007 sample. 


Econometric tests of the hypothesis that January 1983 and December 1998 constitute two breakpoints provide mixed evidence. 
Running Chow breakpoint tests and CUSUM tests on the individual equations of the vector autoregression Eq. (17), we find that 
the null hypothesis of no structural break is rejected for several equations (see Table 7) at the 95% confidence level. The rejections 
are more frequent with the Chow test than with the CUSUM test, probably indicating that, although coefficients change after the 
breakpoints, the predictive accuracy of the regressions does not deteriorate much by extending the sample period. However, if we 
run a LR test of the joint hypothesis of stability of all the equations in the VAR (we use the LR statistic adjusted for small-sample 
bias proposed by Sims, 1980 and reported by Hamilton, 1994), we are not able to reject the hypothesis of no structural change (the 
p-values are 67 and 73% for the pre-1983 and post-1998 periods respectively). Furthermore, the majority of the empirical findings 
presented in the previous sections remain broadly unaffected by the enlargement of the sample period. We now briefly report the 
differences with respect to the base sample (1983-1998). 

Estimating the model with the 1974-1998 sample, we observe the following features: pricing accuracy only slightly 
deteriorates, with the average pricing error remaining well below 10 basis points; the ability of model-implied risk premia to 
predict excess returns on bonds is comparable to that found in the 1983-1998 sample; the R? found in the regressions of 
realized on model-implied expected excess returns is on average higher than 30% and the p-values found in the tests for 
unbiasedness are even greater than those found with the baseline sample; we still find evidence of countercyclicality in risk 
premia, although less pronounced; variance decompositions of the exchange rate indicate that risk premia still explain an 
important fraction of exchange rate variability, although the explanatory power increases at short and medium forecasting 
horizons and decreases at longer horizons; impulse response functions reveal that the delayed overshooting phenomenon 
remains, both for short term interest rates and for risk premia, although it is somewhat weaker for short-term interest rates; 
deviations from uncovered interest rate parity are still driven by positive return differentials both in the short and in the long 
segment of the bond market. 

Estimating the model with the 1983-2007 sample, we find that: pricing accuracy is virtually unaffected; the ability of model- 
implied risk premia to predict excess returns on bonds is comparable to that found in the 1983-1998 sample for German bonds 
while it slightly deteriorates for US bonds (both from the perspective of a German and of a US investor); the R? found in the 
regressions of realized on model-implied expected excess returns is however higher than 25% also for US bonds; the tests for 
unbiasedness still yield very satisfactory p-values; countercyclicality in risk premia disappears; performing variance 
decompositions of the exchange rate, we find that risk premia explain an even greater fraction of exchange variability, due to 
the fact that also German risk premia now play a non-negligible role (remember that in the 1983-1998 sample only US risk premia 
seemed to matter); estimates of impulse response functions reveal that the delayed overshooting phenomenon remains an 
important phenomenon, both for short term interest rates and for risk premia; deviations from uncovered interest rate parity are 
still driven by positive return differentials both in the short and in the long segment of the bond market. 

Overall, extending the sample backwards and forward we do not find striking differences in the results, except for the fact that 
adding recent data takes the evidence of countercyclicality in risk premia away. The findings that a portion of the variability of the 
exchange rate are explained by risk premia and that delayed overshooting happens also in response to changes in risk premia, two 
of the main empirical findings of the paper, are robust across samples. 


4.9. Other robustness checks 


Apart from experimenting with different samples, as documented in the previous subsection, we performed a number of other 
robustness checks both on the data set and on the specification of the model. 

We replaced HP-filtered industrial production with the 12-month growth rate of industrial production and found no significant 
differences. We also computed a monthly series for real GDP by linearly interpolating the quarterly observations and we applied to 
it the same HP filter used for industrial production: although the series based on industrial production is more volatile, it is highly 
correlated with that based on GDP and assigns approximately the same dates to peaks and troughs of economic activity. 
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We added more unobservable variables to the model, but found no appreciable improvement in pricing accuracy or differences 
in the dynamics of the model. We added up to four lags of the macroeconomic variables and the exchange rate, and we found no 
significant changes in variance decompositions and impulse response functions; however, we noted that adding more lags caused a 
fast growth in the condition number of the data matrix and we thus chose the most parsimonious specification (only one lag) to 
avoid collinearity problems. 

We tried different Choleski orderings of the variables for performing variance decompositions: although numerical results 
obviously change, none of the main features reported in Subsection 4.5 are substantially altered. 

We estimated a vector autoregression like Eq. (17), but excluding the two variables related to bond risk premia, and we observed that 
the proportion of variance of the exchange rate left unexplained increased by more than 10% at all time horizons. This check should 
ensure the robustness of the finding that a significant portion of exchange rate variability is driven by bond risk premia. 

We ran Dai and Singleton's (2002) LPY tests to assess the ability of the model to reproduce both the historical and the risk- 
neutral dynamics of yields and we found no evidence of mis-specification. 

Finally, we checked for structural stability during the 1983-1998 sample and we did not come across any obvious structural change. 


5. Conclusions 


We have proposed a no-arbitrage term structure model that allows to contemporaneously price bonds in two different countries, 
taking into account the dynamics of the exchange rate and of observable macroeconomic variables such as inflation and the output gap. 
In the model, the domestic pricing kernel and the exchange rate are specified exogenously and the foreign pricing kernel is derived 
endogenously. Thanks to this specification strategy, the model affords a considerable analytical tractability and allows to analyse 
bidirectional linkages between bond prices and the exchange rate. Estimating the model with US and German data over various sample 
periods, we have found that a significant portion of the variability of the exchange rate is accounted for by time-varying bond risk premia 
(more than 20% over long time horizons). We have analyzed the dynamics of the estimated model to seek explanations for this finding. 
Impulse response analysis reveals that a currency tends to persistently appreciate when risk premia on long-term bonds denominated 
in that currency rise, that is when investors expect large capital gains on long-term bonds denominated in that currency. The delayed 
overshooting phenomenon found by many previous studies with reference to short-term policy rates seems to extend also to expected 
returns on long-term bonds. Furthermore, differences in bond risk premia between countries drive deviations from uncovered interest 
rate parity: the higher is the difference between bond risk premia in two countries, the more profitable is a carry-trade strategy based on 
such difference. After controlling for macroeconomic variables, bond risk premia in one country have fairly limited influence on those of 
another country. Finally, exogenous shocks to the exchange rate have a negligible impact on the yield curves. 

While the econometric specification and the estimation techniques we have proposed allow to efficiently tackle the high 
dimensionality and the computational complexity of a two-country no-arbitrage model, we have constrained ourselves to a 
Gaussian and linear setting. Since there is abundant evidence of non-linearities both in exchange rate dynamics (e.g.: Mark & Moh, 
2005; Nirei & Sushko, 2011) and in interest rate pricing (Ahn, Dittmar, & Gallant, 2002), it could be interesting to explore how the 
model could be improved to take this empirical evidence into account. For example, one could relax the linearity assumption and 
allow for quadratic terms in the pricing equations, along the lines suggested by Ahn et al. (2002) and Diez de los Rios (2009). We 
leave these kind of extensions for future research. 


Appendix A 


A.1. Pricing formulae for foreign bonds 


We derive the pricing formulae for foreign bonds recursively. The price of a bond expiring in one period is: 
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Hence, 
zi = —In(q) 
Tis Mot T_ f= T (22) 
= (ofa sf ser) + (ofp + f')X 
so that: 
1 = 
C, =a-f'p- 5f 22S D, =b-p'f +f (23) 
For foreign bonds expiring in n periods, we have: 
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A.2. Measurability of the currency depreciation rate 
The standard argument used to derive the currency depreciation rate endogenously, as a function of the domestic and foreign 


pricing kernels, proceeds as follows (e.g.: Backus et al., 2001). Let be the payoff at time t+ 1 to an asset denominated in the 
foreign currency. Its price p(§), can be derived either using the domestic pricing measure Q: 


(25) 
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or the foreign pricing measure Q’: 


(26) 


t 


where r{ =z} is the foreign short rate. Assuming that the two above pricing relations hold for any asset €, which is P-mesurable and 
square integrable with respect to P, the equality 
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thanks to a basic fact from Hilbert space theory (see e.g. Luenberger, 1969, page 48, lemma 2). Hence, the currency depreciation 
rate is equal to the ratio between the foreign and the domestic pricing kernels. 

There is a subtlety involved in the above derivation. In order for the expected value in (25) to have any meaning at all, the 
currency depreciation rate S,,4/S, must be an F-measurable random variable. This is by no means an innocuous assumption. 
Recent macro-finance studies (e.g.: Chabi-Yo & Yang, 2007; Dong, 2006) define the probability space (Q, F, Fe P) implicitly by 


defining the law of motion of ;: F+ is the filtration generated by X+, F = o( U, F) ,and Pis the product measure derived from the 
= 


transition densities of X; hence, assuming that the currency depreciation rate be F-measurable is tantamount to saying that the randomness 
in S;is completely explained by the randomness in X;. While this is obviously true if X; explicitly includes the currency depreciation rate, it is 
hardly a realistic assumption otherwise, as shown by the empirical evidence provided by, among others, Inci and Lu (2004). 


A.3. Further identifying restrictions: real Schur decomposition of the autoregressive matrix 


Given the following parametrization of the model: 
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the set of identifying restrictions proposed by Pericoli and Taboga (2008) is: 


1. $° is lower triangular; 
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2. >"°=0; 
CEP ae 
4. b">0; 
5. Xp=0. 


As pointed out by Hamilton and Wu (2010), these restrictions are not sufficient for identification: in fact, premultiplying X¥ by 
an orthonormal matrix H, we can obtain an infinite number of equivalent representations of the system that satisfy the same 
parameter restrictions: 
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where the new vector of unobservable variables is: 
Xe SX, 


The innovations to the unobservable variables are: 


e™" — He" ~N(0, Ipu) 
ne" = He ~N(0, le) 


because 


Vare "] Hla H" = HH" = Ip 


and by the same token Var [nf “| = Ipu. 

The new parameters are: pt u— Hp", p™ OU UH. u” “= Hu, pr 40 — Hpo, pt uu — HoH", pie = pH", m = Ho", 
p” = Hp”, paw = Hp""H". 

To avoid this lack of identification, we can use a well-known theorem on real Schur decompositions. 

Define the real Schur property as follows: 


Definition 1. A kxk matrix T satisfies the real Schur property if it is a block upper triangular matrix whose entries Tj satisfy the 
following restrictions: 


1. the diagonal entries of T are equal to the real parts of the roots of a characteristic equation of order k, i.e.: 
T; = re(A;) io eee 
where A; (i= 1, ..., k) are the solutions of a characteristic equation of order k: 
Faa beh +. iG =O 
with real coefficients Ck—1, Ck— 2, ..., Co, and 
fi > To ean Ty 
2. the entries below and above the principal diagonal satisfy: 
Tii-1Ti-11 = ref) - re(Ni-1) ~N Nai = 2,...,k 


and 


Tij = 0 if re(Aj)- re(Ay4) Aj Ai- = 0 
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3. the entries below the diagonal blocks are zero: 
Ty =0 if i>j+1 


A standard result in matrix algebra, known as the real Schur decomposition, states that any real kxk matrix A can be 
decomposed as: 


A=H'TH 


where H is a kx k orthonormal matrix and T is a kx k matrix satisfying the real Schur property (see e.g. Bernstein, 2005 - Theorem 
5.4.1 and Corollary 5.4.2). 

The real Schur decomposition can be used to avoid the identification problem mentioned above. This can be done by imposing 
that the rotation matrix be equal to the matrix H in the real Schur decomposition of p™". In other words, when p™ does not already 
satisfy the real Schur property (in which case the only possibility is H=/), we pin down H by requiring that the rotation 


pum Z HoH" 


transforms p“” into a matrix p™ “ satisfying the real Schur property. The requirement of increasing diagonal entries makes sure 
that no further rotations (different from H=J) are possible when p™ satisfies the real Schur property. Therefore, the lack of 
identification can be resolved by imposing the restriction that p satisfies the real Schur property. 

Also note that the constraints imposed by the real Schur form are such that any 2 x 2 diagonal block of p"“ can be written as a 
function of 3 parameters, because the entry in its lower left corner is automatically determined by the other 3. Therefore, there are 
(ku +1) xk,/2 free parameters in p““, which matches the parameter count for exact identification done by Hamilton and Wu 
(2010). The second potential identification problem pointed out by Hamilton and Wu (2010), which might come from the 
presence of unit roots, is avoided by the identification condition X§=0. 


A.4. Expected excess holding period returns 


Remember that the expected holding-period excess return on a domestic bond (from the perspective of a domestic investor) 
is: 


eS Er [In(p X'h) —In(pr) | i in(pt) 


Given that bond prices are exponential affine functions of the state variables, we obtain, by substituting the pricing formulae of 
Section 2: 
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—(n—h) (An-n + By-nEt [X + al ) H n(A, H B,X;) h(A; t BX) 


Since 


E [X +h] =H+ PX ha 


and 


by recursive substitution one obtains: 
Eş [Xe + n] = Yh + p"X, 
where 
v, =u + pu +... tpn 
Hence, 


„h 
e; = Jah e KanXı 
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where 


Jan = = (An-n + Bla) + nAn —hA, 


Kin = —(n—h)B; np" + nB} —hB; 
Analogous algebra yields the expressions for the holding-period return on a foreign bond (from the perspective of a foreign 
investor): 


i = Qnn + WanXe 


where 
Qnn = —(n—h) (Ca-n + Da—hYh) + nc, —hC, 
ah = —(N—h)D;_,p" + nD; —hD; 
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